Nonequilibrium work statistics of an Aharonov-Bohm flux 

Juyeon Yi 

Department of Physics, Pusan National University, Busan 609-735, Korea 

Peter Talkner, Michele Campisi 

Institute of Physics, University of Augsburg, 
Universitdtsstrasse 1, D-86135, Augsburg, Germany 
(Dated: June 28, 2011) 

Abstract 

We investigate the statistics of work performed on a noninteracting electron gas confined into 
a ring as a threaded magnetic field is turned on. For an electron gas initially prepared in a 
grand canonical state it is demonstrated that the Jarzynski equality continues to hold in this 
case, with the free energy replaced by the grand potential. The work distribution displays a 
marked dependence on the temperature. While in the classical (high temperature) regime, the 
work probability density function follows a Gaussian distribution and the free energy difference 
entering the Jarzynski equality is null, the free energy difference is finite in the quantum regime, 
and the work probability distribution function becomes multimodal. We point out the dependence 
of the work statistics on the number of electrons composing the system. 

PACS numbers: 05.30.-d, 05.70.Ln, 05.40.-a 
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I. INTRODUCTION 



One of the most fascinating electro-magnetic field effects is the modulation of quan- 
tum interference in multiply connected spatial regions due to electromagnetic fields. The 
Aharonov-Bohm (AB) effect is a well-known example where a localized magnetic field intro- 
duces the phase shift of a particle wavefunction and results in an interference pattern gov- 
erned by the AB flux ^ab = / A ■ cil A similar effect, called the Aharonov-Casher (AC) 
effect occurs, when a neutral particle with a magnetic moment /x moves in an electric field 
E and acquires a phase shift amounting to the flux ^ac = §^^d\- {jl/c j24|. A dual effect 
was also pointed out for a neutral particle carrying an electric dipole moment moving in a 
magnetic field of the appropriate configuration |5| . In the mentioned examples, electromag- 
netic fields influence the wave function and also the energy spectrum of a particle moving 
in a multiply connected spatial region but do not exert any classical force. 

Recently, scientists working in the field of non-equilibrium thermodynamics have drawn 
the attention to the fact that the work done by external forces on a driven system may 
be usefully employed to characterize their response properties. Jarzynski introduced the 
celebrated nonequilibrium work relation that links the free energy difference (AF) to the 
averaged exponentiated negative work {g, 7|: 



(1) 



where w is the work performed on a system by a time- dependent force determined by a 
prescribed protocol and (■ • ■ ) denotes the average over many realizations of the forcing 
experiment. The equality was primarily derived for classical systems, which experiments 



and theories so far mainly refer to js-llOl. Fluctuation theorems in presence of rnagnetic 
fields and other non-conservative forces were studied for classical systems in Ref. On 

n 

the other hand, generalizations of Eq. (1) to quantum mechanical systems were discussed |12h 
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221 ]. for recent reviews see 



24( 1 . In quantum mechanics, the work is obtained by means 
of two energy measurements at the beginning and at the end of a given protocol. In the 
mentioned examples of quantum interference effects, the electromagnetic fields do work on 



a charged particle, or a magnetic or electric dipole in multiply connected domains not only 



caused by the c 



spectrum |25l-l27|. 



assical forces exerted on the particle but also by the shifts of the energy 
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The main purpose of tliis work is to investigate the statistics of the work done by an 
Aharonov-Bohm flux for quantum charged particles moving along a one-dimensional ring 
representing the simplest possible multiply connected domain. We consider not only the 
single-particle case but also many-particle systems by generalizing the fluctuation theorem 
for a grand canonical initial state. Especially we focus on fermionic systems in a ring- 
configuration (see Figj_l). Their equilibrium properties were discussed in terms of persistent 



currents decades ago 



25- 



Recent measurement using nano-cantilevers to detect changes 



in the magnetic field produced by the current has achieved high accuracy and prompted 



a renewed interest in this topic 



29|. It is noteworthy that the focus of our study is laid 



upon the nonequilibrium nature of the system, revealed in the statistics of work done by the 
magnetic flux. By obtaining an analytic expression of the characteristic function of work, 
we examine the quantum and classical nature of the resulting distributions and study their 
dependence on both temperature and particle number. 

The paper is organized as follows: Sec. II is devoted to the introduction of the system of 
interest. In Sec. Ill, we obtain the probability distribution for a single-particle case. We then 
consider a grand canonical initial state of many-particle systems, and present the results in 
Sec. IV. Summary and conclusion are drawn in Sec.V. 
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FIG. 1: (Color online) The left panel gives a schematic picture of a ring with a threading flux, 
and the right panel displays the work distribution consisting of a series of peaks for a single particle 
moving along the ring upon a sudden switch of the dimensionless flux / = ^/^o from / = to 
/ = 0.1. The thick solid line represents the Gaussian approximation of the work distribution, 
Eq. dSD, for the same parameter values. 
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II. THE SYSTEM 



We consider spinless fermions moving along an infinitely thin ring of radius R in 



28 



presence of a magnetic flux. The corresponding Hamiltonian reads 

where 6'^ is the angular coordinate of the £th particle, and / = $/$o is the total flux 
threading the ring, $, in units of the flux quantum, $0 = hcje. The single particle energy 



eigenvalues are given by 



28 



E,(/) = Eo(A;-/)^ (3) 

where -Eq = f? jlvnl^ characterizes the energy-level spacing, and the integer k denotes the 
angular momentum quantum number. 

For later consideration of a many-particle system, let us introduce the second quantized 
form of the Hamiltonian: 

= Y,Ek{f)cW = 5^i?fc(/)A4, (4) 

k k 

where (cfc) is the creation (annihilation) operator of an electron in the kih. angular mo- 
mentum (or energy) eigenstate, and the number operator Hk = c^C/t measures the particle 
number in the /c-th state. 

We will study the probability distribution function P{w) of the work w that is performed 
on the electrons in the time span [0,r], by a magnetic flux f{t) that varies in time. As a 
consequence, the Hamiltonian of the system becomes time-dependent. It will be denoted by 
'H(t) = 'Hf(t)- Here we restrict ourselves to the case of a sudden switch of the magnetic flux 
immediately after the time t = with 

{0 for t < 
- . (5) 
/ for t > 

Given this protocol, we will first calculate the work characteristic function (i.e., the Fourier 



transform of P{w)), which for a canonical initial state is given by the formula 16l-|l8l |21| : 

00 

Giu) = / dwe''''"Piw) = (e»««We-^"^(o))p^, (6) 
J 00 
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and then obtain P{w) by inverse Fourier transformation. Here {X)p^ = TrXpc and 
Pc = e'^'^^^y Zq with the normahzation Zq being the canonical partition function. Fur- 
ther, 'HH{t) denotes the Hamiltonian operator in the Heisenberg representation. Note that 
since ['H(t), 'H(s)] = for any t and s in the time span, the Hamiltonians in the Schrodinger 
and Heisenberg picture coincide, i.e. l-iHij) = 'H(r). 



III. SINGLE PARTICLE CASE 



In the case of a single particle system, we obtain from Eqs. ([2]) and ([5]), 

G{u) = e'"^°^' Yl e2^'=/^o"e-^^«'=7^o- (7) 

k 

We rescale all variables with the dimension of an energy by the natural energy unit Eq as 
u = Equ, w = w/Eq, and /3 = (3Eo (for notational simplicity, we drop the tilde in the 
following). Note that at high temperatures the system enters the classical regime. In this 
regime we can discard the discretness of k and replace the summation by an integration, 
namely, (R/h) J^^dp with the momentum defined hy p = hk/R. We thus obtain 

G{u) ^ Gc{u) = e*-^^""-'^^*'^/^, which leads to a Gaussian distribution of work reading 



Note that for u = if], Gc{iP) = {e~^^) = 1. This conversion from discrete summation 
to integration is accurate only at sufficiently high temperatures (/3 ^ 1), or for a ring of 
sufficiently large radius. When we fully count the level discretness, the work distribution is 
expected to be a series of peaks, and the normal distribution would provide its envelop. This 
can be confirmed by evaluating Pc{w) directly from Eq. ((71) via an inverse Fourier transform: 

Piw) = J2^k6iw-2kf-f) (9) 



k 



where Wk = e ^'^'^ / e ^^"^ is the weight of the fcth peak. The right panel of Fig. 1 displays 



the resulting distribution 



30|. It is then a due course to investigate many-particle cases. 



where the effect of finite particle number comes into question. 
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IV. GRAND CANONICAL INITIAL STATE 



In order to deal with many-particle systems, we extend the characteristic function of work 
to initial grand canonical states. Although not needed here, we allow for possible changes 
of particle numbers. This will lead to a generalization of the Jarzynski work theorem to 
grand canonical initial states. Only later we will specify to the case of strict particle number 
conservation. 

In the grand canonical ensemble energy and particle number are fluctuating quantities. 
In order to determine their changes effected by a protocol a simultaneous measurement of 
energy and particle number must be performed at the beginning and at the end of the pro- 



tocol, see Ref. [3l| where the joint statistics of changes of two observables are discussed. 
Joint measurements necessitate that the particle number operator M commutes with both 
the initial and final Hamiltonian, i.e., [A/", 'H(r)] = [A/", 'H(O)] = 0. Then joint eigenfunc- 
tions |\l/jy,Ar(0)) and \^u,n{'t)) exist with corresponding pairs of eigenvalues E^, N and E^, 
iV, satisfying H{t)\^ ^^n{t)) = Ep,^(r)|^'p,^(r)) and M\^p,n{t)) = N\^^^f^{T)), as well 
as 7{(0)|^^,jv(0)) = E^(0)|^^,jv(0)) and N\^^^n{^)) = Ar|^^,7v(0)). The joint probability 
density function V{w,n) to observe the work w and particle number change n in a single 
realization of the protocol is given by 

V{w, n) = J2Y. -^'^(^) + EumS^,N-N (10) 

i^,!^ N,N 

X P{E,{t),N\EM,N)P^"''\EM,N), 

where P^^'^^ = / Qq is the joint probability of finding the energy -E'^(O) and particle 

number in the initial grand-canonical state; further, Qq = Xle^ ^ ^-l^i^^i^)-!^'^) that is the 
grand-canonical partition function. The conditional probability P{Ep{T),N\Ey{{]),N) for 
finding the energy Ep{t) and particle number N at the end of the protocol once it was £'^(0) 
and at the beginning is determined by the overlap between the final state and the time 
evolved initial state: 

P{E,{t),N\E,{0),N) = |(vl/,(r)|t/(r,0)|*,(0))|^ (11) 

where U{t, 0) is the unitary time evolution operator solving the Schrodinger equation 
ihdU{t, 0) /dt = n{t)U{t, 0) with f/(0, 0) = 1. The Fourier transform of the joint probability 
( 1T0|) then yields the characteristic function that can be cast into the form of a two-time 
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correlation function 3l|, i.e. 



oo „ 



(12) 



Ju'HH{T)+ivNu{T) p-iu'H{0)-ivN{ld) 



Pa- 



Putting w = if3 and v = —i(3fi one obtains a generalized Jarzynski equality for the grand- 
canonical initial state reading 

Qr 



(13) 



with the grand potential difference AQ = fi(r) — Q{0) where Q(t) = — /3 ^lnQ(t). Similar 
considerations were made for classical system s [32 1 and also for composed quantum systems 



with number exchanges between subsystems [3^ 



35|. 



V. MANY-PARTICLE CASE 



We analyze the work statistics of a many electron system undergoing a sudden switch 
of the magnetic flux by means of the generalized Eq. f|T3|) . Since in our case the particle 
number is a constant of motion, A/i^(r) = A/'(0), the characteristic function, Eq. ( !T2|) . is 
independent of v; therefore we simply write it as g{u, v) = Q{u). Due to the sudden switch 
of the magnetic flux 'Hi^(r) = 'H(r) = J2k Ek{,f)-^k as given by Eq. (3). Moreover 7/(0) and 
'H(r) commute with each other, and we can then write 

Q{u) = (e»S'^-^^(/)^'^-),, = nil - (A4) + (A/'fc)e»^'^-(^)], (14) 

k 

where we used the property JV^ = Mk of fermionic number operators. Here Afc(/) = 
Ek{f) — Ek{0) and (A/1-) = [1 + e'^'^'^*^*-'^''"'''']"^ for the fermionic particles. The chemical 
potential /i should be determined to satisfy 

k 

where (A/") is the average number of particles, which will be denoted as hereafter. 

Figure 2 shows the work distributions for the flux / = 0.1, at different temperatures and 
particle-numbers. As shown in Fig. 2(a) at low temperatures (large /3), the distribution is 
narrow and centered at the difference between the ground-state energies in presence and 
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FIG. 2: (Color online) The probability distribution 'P{w) of the collective work performed on N 
spinless fermions by a sudden switch of the magnetic flux from / = to / = 0.1. In panel (a) 
results are displayed for the average number of particles = 25 at three different temperatures. 
The solid line represents the high temperature Gaussian approximation Eq. for /3 = 0.02. Panel 
(b) displays the distributions for two different average particle numbers at inverse temperature 
P = 0.1. Note that for the larger N value the side peaks have a larger distance from the central 
peak. 



absence of the magnetic flux which is given by Wc = 'Ylk&kgXf'^ ^ 2A;/). Here the summation 
runs over the /c-values given by kgg = 0, ±1, • ■ ■ , ±(A^ — l)/2 for odd A^. Due to the pairwise 
cancelation of positive and negative /c-values in Wc, the term linear in / vanishes, and hence 
Wc = Nf^. For the used parameters = 25 and / = 0.1, Wc = 0.25 which indeed coincides 
with the central peak positions in Fig. 2(a). At higher temperatures (see /3 = 0.1 in Fig. 2a), 
excited states ofk = ±(A^+l)/2 come into play, which lead to side peaks located at Wg = Wc=t 
(A^+ 1)/. The number dependence of the side peak positions can be seen in Fig. 2(b): With 
decreasing particle numbers the distances between the central peak and the side peaks shrink. 
At high-temperatures, many excited levels contribute to the work fluctuation with almost 
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equal weights. This leads to the seemingly continuous and flat distribution as displayed for 
P = 0.02 in Fig. 2(a). In this case, in fact, particles follow the Maxwell-Boltzmann statistics, 
and the energy spectrum can be conceived as a continuum. Then, the characteristic function 
is approximately given by the products of single particle contributions, i.e., by Qc{u) ~ 
G^{u). This gives Vc{w) ^ ^ I3/{A-itN p)e-^^'"-^f^^^/^^^f^^ which is plotted as solid line for 
P = 0.02 in Fig. 2(a). We present the temperature dependence of the probability distribution 
in Fig. 3(a). It displays the change from a narrow unimodal distribution at low temperatures 
through a multiply peaked distribution at intermediate temperatures to a broad Gaussian 
distribution at high temperatures. From the characteristic function of work one obtains the 
variance and all nth-order cummulants, C„ via the formula, C„ = (— i)"'<9" ln^(M)|„=o; where 
9" denotes the nth derivative with respect to u. As shown in Fig. 3(b) the variance and the 
third order cummulant for A^ = 25 rapidly decrease to zero with decreasing temperature. 
The inset shows the exponential temperature dependence of C3. We note that also the 
variance cr^ = C2 decays exponentially with temperature, although we do not shown it here. 
On the other hand, for A^ = 26 the variance saturates to a finite value, while the third 
moment exhibits an exponential decay, similarly to, but far faster than for the case when 
A^ = 25. 

The explicit form of the variance for the system is given by 

<r' = J2^{i-{mwAi{f)] (16) 

k 

which indicates that the number fluctuations (1 — (A/fc))(A/fc) determine the variance of the 
work. The total average number A^ of particles in the initial equilibrium state in absence of 
a magnetic flux is given by Eq. ( fT5l) . At low temperatures it increases in a stepwise fashion 
with varying fi and forms plateaux of height A^ = 2n + 1 with steps close to = n^ (see 
Fig. 4(a)). This behavior is a direct consequence of the degeneracy of states with angular 
momentum ±k. Due to the jumps, a system with an even average number of particles 
has pronounced work fluctuations that persist with decreasing temperature whereas for an 
odd number the variance of work vanishes with decreasing temperature. We present the 
number dependence of the standard deviation in panel (b) of Fig.4. At the low temperature 
(/3 = 5.0), a vanishes at odd A^ (except A^ = 1), whereas it has peaks at even A^'s. As the 
temperature increases, the dips and peaks at small and intermediate average numbers N 
merge into a smooth and increasing curve but remain visible at sufficiently large values of 
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FIG. 3: (Color online) (a) The probability distribution of work w for = 15 under a sudden 
magnetization flux switch from / = to / = 0.1 is displayed as a function of w and the inverse 
temperature /3 by means of different gray- values as specified in the left upper part of the panel. At 
low temperatures the distribution is unimodal and develops side peaks with increasing temperature, 
(b) The temperature dependence of the standard deviation of work a is compared for N = 25 and 
N = 26. In the case of the even average particle number the standard deviation saturates at a 
finite value with decreasing temperature while it vanishes for the odd average particle number. 
The inset shows the natural logarithm of the third cummulant (C3) versus /3. 

N. 

VI. SUMMARY AND CONCLUDING REMARKS 

In summary, we investigated the work distribution of many non-interacting fermionic 
particles driven by an AB flux in a non-simply connected geometry. In the single-particle 



10 




-5 5 10 15 20 25 30 



— I — I — I — I — I — I — I — I — I — I — I — I — 1 — I — 1 — I — I — I — I — I — 1 — I — I — I — I — 




N 

FIG. 4: (Color online) (a) The average number of particles as a function of the chemical potential 
for various temperatures. The stepwise increase of N can be seen at low temperatures (see the 
curve for /3 = 5.0 and the inset for /3 = 20). (b) The number dependence of the standard deviation 
of the work at / = 0.1. At the low temperatures care displays narrow dips at odd N and broad 
peaks at even N. At the lowest temperature /? = 20 shown in the inset the standard deviation is 
vanishingly small up to = 1 and then displays a maximum at A^ = 2, see the inset. At high 
temperatures (/3 = 0.3), those structures are washed out at small and intermediate A^ and become 
visible only for sufficiently large A^. 



11 



case the work distribution at high-temperatures, namely, in the classical regime, is given by 
a Gaussian distribution, yielding {e~'^^)=l, indeed confirming that the quantum flux leaves 
the free energy unchanged. By contrast, the distribution in the quantum regime is found to 
be multimodal caused by particle excitations. In particular, in order to deal with a many- 
particle system, we have generalized the expression for the characteristic function of work to 
quantum systems that initially are in a grand canonical state. We proved that the difference 
of the grand potentials of a hypothetical grand canonical equilibrium system with the initial 
temperature and chemical potential at the at the flnal parameter values and of the actual 
initial system enters a generalized Jarzynski equality. Although an energy measurement 
is an experimentally challenging task, theoretical examination of work in quantum many- 
particle systems per se is worthwhile for the fundamental understanding of nonequilibrium 
characteristics. 
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